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ABSTRACT 

Using our non-local time-dependent theory of convection, the linear non-adiabatic 
oscillations of 10 evolutionary model series with masses of I-SMq are calculated. 
The results show that there is a red giant instability strip in the lower temperature 
side of the Hertzsprung-Russel (HR) diagram which goes along the sequences of the 
red giant branch (RGB) and the asymptotic giant branch (AGB). For red giants of 
lower luminosities, pulsation instability are found at high order overtones, the lower 
order modes from the fundamental to the second overtone are stable. Towards higher 
luminosity and lower effective temperature, instability moves to lower order modes, 
and the amplitude growth rate of oscillations also grows. At the high luminosity end of 
the strip, the fundamental and the first overtone become unstable, while all the modes 
above the 4th order become stable. The excitation mechanism have been sdudied in 
detail. It is found that turbulent pressure plays a key role for excitating of red variables. 
The frozen convecttin approximation is unavailaable for the low temperature stars with 
extended convective envelopes. In any case, this approximation can explain neither the 
red edge of the Cepheid instability strip, nor the blue edge of the pulsating red giant 
instability strip. An analytic expression of a pulsation constant as a function of stellar 
mass, luminosity and effective temperature is presented from this work. 
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1 INTRODUCTION 

The Mira type variables are among the most luminous stars 
in a star cluster, and they possess very large light ampli- 
tudes therefore they are the first variable stars ever discov- 
ered. The red variables, however, are not yet understood 
well. According to Wood (2000) "red giant stars are proba- 
bly the least understood of all variable stars" . It was suspi- 
cioned long time ago that all red stars are variables. Eggen 
(1969) pointed out, for instance, that all old disc popula- 
tion red giants with (R-I)o > 1-0 are variables. A few year 
later, he shifted the blue edge of that instability area to (R- 
I)q « 0.9 (Eggen 1973). The visual amplitudes of red vari- 
ables have vary large ranges, from 0.01 to 10 magnitudes, 
and the typical time scales (periods) goes from several to a 
few hundred days. Classification of these red variables based 
on period, luminosity and population data were made by 
Eggen (1972a,b, 1973a,b, 1975, 1977). By photometric and 
spectroscopic observations of 187 G & K type giants, Henry 
et al. (2000) showed that all K5-M0 type giants are vari- 
ables, and they further claimed that the variations of red 
giants with spectral type later than K2 are resulted from 
pulsation. Towards lower effective temperature and higher 
luminosity, the light amplitude tends to increase. In fact, 
these stars are located either in the upper part of RGB or 



on AGB. It is identified by Percy et al. (1998, 2001, 2003) 
that these stars are actually pulsating from the fundamen- 
tal up to the 3rd overtone. It is also found that there are 
variables with rather short periods from several hours to ~2 
days and small amplitudes from several to 20 percents of a 
magnitude on the lower RGB of globular clusters (Yao & 
Qin 1993, Yao, Zhang & Qin 1993). By analyzing the 38.5 
day photometric data acquired by HST for 47 Tuc, Edmods 
& GUiland (1996) discovered a mount of small amplitude K 
type variables on RGB, whose periods go from 2 to 4 days 
and the amplitudes only measure 5-15mmag. On the V-(B- 
V) color magnitude diagram (CMD) of 47 Tuc, one group of 
these variables gathers near the red edge of the horizontal 
branch (HB), while another stays on the RGB and the bot- 
tom of the AGB about 1 magnitude brighter than the HB 
of the cluster. 

The observational studes of red giant variables gains are 
advanced by the variable star survey programs and gravita- 
tional lensing survey projects started in recent years, such 
as the Long-Term Photometry of Variables (Jorissen et al. 
1997), the All Sky Automated Survey (Pojmanski 2002) and 
the Northern Sky Survey (Wozniak et al. 2004) etc., es- 
pecially, MACHO (Alcock et al. 2000, Wood et al. 1999, 
Schultheis 2004), OGLE (Kiss & Bedding 2003, 2004, Wrag 
et al. 2004, Soszynsky et al. 2004), EROS (Cioni et al. 2001) 
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and MOA (Bond et al. 2001, Noda et al. 2002) have accu- 
mulated a huge time sequenced stellar photometric data set 
that is suitaable for searching of variable stars. Wood (2000) 
found 5 distinct period-luminosity sequences (A, B, C, D & 
E, as he defined) in the K-log P diagram based on his analy- 
sis of the MACHO near infrared J & K broadband photome- 
try data in the half square degree field of the LMC bar from. 
In that, the sequence C was the Mira variables identified as 
the fundamental radial mode. The stars on the sequence 
B wooud be known as semi-regular (SR) variables and A 
are of very low amplitude as small amplitude red variables 
(SARVs). The sequence B and A correspond to pulsation in 
first-third overtines (Wood, 2000). Such result is later con- 
firmed independently by OGLE observations (Kiss & Bed- 
ding 2003, 2004). It is further pointed out by Soszyriski et 
al. (2004) that the OGLE small amplitude variable red gi- 
ants (OSARGs) and the long period variables (LPVs) are 
two distinct types of pulsating red variables. 

The excitation mechanism of the pulsating red variables 
are far from fully understood, the principal barrier is that 
there is not yet a perfect time-dependent convection the- 
ory that can handle the coupling between convection and 
oscillations. Most of the studies either simply neglected the 
coupling between convection and oscillations, such as Ostlie 
& Cox (1986) and Gong et al. (1995), or over-simplified it, 
such as Cox & Ostlie (1991), Fox & Wood (1982). The most 
careful treatment of the coupling between convection and 
oscillations comes from Balmforth (1992) and Xiong, Deng 
& Cheng (1998), with the former study applied Cough's 
(1977) non-local time-dependent mixing- length theory of 
convection, while the later one adopted Xiong's (1989a) 
time-dependent statistical theory of convection. In our case, 
Xiong, Deng & Cheng (1998) studied the oscillations of Mira 
variables. The results show that for the luminous red giants, 
the oscillations tend to happen in the fundamental and low 
overtones whose amplitude growth rate is very large, i.e. 
the amplitude gets doubled within a few, or sometimes even 
within a single period. All the overtones above 4th order 
are stable. For intermediate and low luminosity red giants, 
however, the fundamental mode and the 1st overtone are sta- 
ble while being unstable at modes above the 2nd overtone 
(Xiong & Deng 2001a). These SARVs have rather small am- 
plitude growth rates, being almost neutrally stable, there- 
fore special care is needed to handle convection in this situ- 
ation. The goal of the present work is to deliver a thorough 
study of pulsation stability in stars with intermediate- and 
low-luminosity on the RGB and AGB. The second section 
describes the refinement and specification of the convection 
theory for the current work, and the third section presents 
the results of linear non-adiabatic computations. The exci- 
tation mechanism for the red variables is discussed in section 
4. The conclusions of this work and discussions are given in 
the last section. 



2 THE FUNDAMENTAL EQUATIONS 

Stars near the Hayashi limit (RGB and AGB) all have very 
deep convective envelope. Instead of the K-mcchanism op- 
erating in the hot and warm stars, the coupling between 
convection and oscillations will be the primaxy mechanism 
of excitation and damping for these cold stars. The main 



difficulty prevent us from understanding the oscillations of 
such stars has been the lack of a perfect convection theory. 
Being distinct from the pulsational properties in the even 
cooler and brighter long period variables, the pulsation am- 
plitudes of the intermediate and low luminosity red giants 
should be very small if any unstable p-mode oscillations can 
be excited. This fact actually requires very careful treat- 
ment of the coupling between convection and oscillations, 
therefore makes the problem especially challenging. 

Convection transfers energy and momentum inside 
stars, and causes mixing of matters. A star, as a thermo- 
dynamic system, can also be regarded as a heating (or cool- 
ing) machine when considering its thermodynamic behavior. 
Based on its nature and mechanism of working, the coupling 
between convection and oscillations can be divided into two 
types as the following: 

(i) The thermodynamic coupling between convection and 
oscillations, whose mechanism is to adjust the temperature 
and pressure of stellar material through its energy trans- 
fer function, therefore affects the pulsational stability of the 
star. Its general efi^ect is damping. The thermodynamic cou- 
pling is the major reason for the existence of the red edge 
of the Cepheid instability strip (including RR Lyrae and S 
Scuti) in the HRD (Xiong, Cheng & Deng 1998, Xiong & 
Deng 2001a). 

(ii) The dynamic coupling between convection and oscil- 
lations, whose mechanism is to trigger momentum transfer 
within a star, i.e. to affect the stability through turbulent 
pressure and turbulent viscosity. Turbulent viscosity is al- 
ways a damping, while turbulent pressure is generally a ex- 
citation for oscillations. Such coupling plays a key role in 
exciting the red pulsating variable stars (Xiong, Deng & 
Cheng 1998). 

The relative contributions of convective energy transfer, 
turbulent pressure and turbulent viscosity to the instability 
of low temperature srars change with temperature, lumi- 
nosity and pulsating mode of the stars, and this fact can 
be used to explains the complex behaviors of pulsation in- 
stability among these stars having extended convective en- 
velopes. It is therefore quite easy to understand such ques- 
tions as follows: Why does the Ceipheid instability strip have 
a red edge? Why do the long period variables exist outside 
the Cepheid instability strip? What makes the luminous red 
variable stars to have higher amplitudes and to pulsate at 
the fundamental and lower order overtones than other stars, 
and what is the reason for the decreasing of amplitudes in 
the red variables with decreasing luminosity and eventually 
becoming pulsationally stable? 

The main challenge in handling the coupling between 
convection and oscillations is the requirement of a time- 
dependent theory of convection. Owing to the huge scale 
of stars, stellar convection is generally fully developed tur- 
bulence which is still not yet understood satisfactorily. The 
main difficukty is due to the nonlinearity and the non- 
locality of hydrodynamics. Compared to the thermodynamic 
coupling between convection and oscillations, the dynam- 
ical coupling is even more difficult to deal with, for the 
anisotropic of stellar convection. The convection theory ap- 
plied in this work is a time-dependent one. Neglecting stellar 
rotation and magnetic field, the full set of equations used in 
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the calculations of stellar convective envelope structure and 
radial oscillations is the following, 
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where r is the radius, Mr is the mass within the sphere 
of radius r, Ur is the radial component of the velocity, p, 
T and P are respectively the regular notations for den- 
sity, temperature and pressure (including radiative), x is 
the opacity, Cp is the specific heat at constant pressure, 
B — — {din p/dlnT) p is the expansion coefficient of gas. 
Our formulism of stellar convection is a statistical theory of 
correlation functions based on hydrodynamics and the the- 
ory of turbulence, x^ a^nd n'-' are respectively the isotropic 
and the anisotropic components of the auto-correlation of 
turbulent velocity w", therefore we have. 



(9) 



Z and V are respectively the auto-correlation of turbu- 
lent temperature fluctuation T' /T and the cross-correlation 
of radial turbulent velocity and temperature, deflned as. 
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In our equations, all the quantities with a prim are the 
turbulent fluctuations, while those with a bar are statistical 
averages of the corresponding quantities. 

Eqs. (H}-® are derived based the equations of hydro- 
dynamics and the theory of turbulence, where Eqs. Q-Q 
are the normal equations for stellar structure, being different 
from the traditional expressions only by having the turbu- 
lent Renold's stress (turbulent pressure and viscosity) in- 
cluded in eq. ^ and having both the turbulent thermal 
convective flux Lc = iivr^pGpTV and the turbulent kinetic 
energy flux Lt = Anr'^ pu'rw[w''^ /2 included in eq. (13} . As 
a result of this consideration, the correlations of turbulent 
convection also show up eqs. ([2]) and (|3}. Eqs. ©-(O are 
a set of dynamical equations of turbulent convection under 
non-local theoretical frame. Owing to the fact that the equa- 
tions of stellar structure include the convective quantities 
and the stellar structural quantities, p, P and T (and Ur as 
well) enter the dynamical equations of turbulent correlations 
of eqs. ([5}-([8ll, it makes these equations coupled with each 
other, therefore eqs. ((T}-(|S} formed a complete set of equa- 
tions for the structure and oscillations of stellar envelope 
at spherical symmetry condition (of course without rotation 
and magnetic fleld, as we assumed at the beginning). In the 
frame work of the local theory, stellar convection is handled 
with a set of algebraic equations, while under non-local ex- 
pressions it is dealt with by a set a differential (or integral) 
equations. In nature, all ffuid motions are in fact non-local, 
the local convection approximation is good only at a limited 
local region (such as the interior zones far away from the 
boundaries of a convective zone). The local theory has no 
chance to be right when considering a whole convective re- 
gion as an entity. This is especially true near the boundaries 
and the overshooting zones where the local approximation 
should be completely lifted, and in this context, one has 
to use a non-local theory. Our non-local theory of stellar 
convection is a statistical theory of turbulent correlations 
derived from hydrodynamics and the theory of turbulence, 
its physical picture is very clear. Eqs. (IS}-(|S)| describe a 
statistical equilibrium of turbulent correlations. Turbulence 
gains energy from buoyant force and the average of motion 
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of fluid, this process happens in the low wave number region 
of the turbulent spectrum. Due to the non-linearity and in- 
teractions of fluid motion, turbulent energy is cascaded to- 
wards higher wave numbers from lower wave numbers, and 
is turned into thermal energy by molecular viscosity and is 
eventually completely dissipated. The right hand side terms 
including the convection parameter ci in eqs. ©-(IS} rep- 
resent here the turbulent viscous dissipation, following the 
turbulence theory the dissipation rate of turbulent kinetic 
energy pe can be expreessed as follows (Hinze 1975), 

e = 2ryefceX^ (12) 

where fee ~ \/3fcei ~ V^/lei is the wave number of the 
energy-containing eddies. It is assumed in this work that Ui 
is proportional to the local pressure scale height, i.e., 

j.2p 
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GMrP 

The third order correlations in eqs. ©-([SJ correspond 
to the non-local convective diffusion. The third order cor- 
relation term in eq. ((5)1, for example, is just the turbelent 
kinetic energy flux. A gradient type of diffusion approxima- 
tion for the third order correlations is adopted in our theory 
of non-local convection. Taking the third order correlation 
for turbulent velocity as a instance, we have. 



Lt = 4-Kr^pu'rWlw''/2 = -Gnr^pxA^^, (14) 

where A is the diffusion length assumed to be propor- 
tional the the local pressure scale height. 
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where C2 is another convective parameter that is linked 
to the non-local turbulent diffusion process. 

Putting eq. (|15|l into eq. (|14|l . Lt can be rewritten as, 

where. 
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Following similar operations, the remaining third order 
correlations in eqs. ((6)|-((8]) can be given as. 
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As described above, in the static station convection 
gains energy from buoyant force, therefore the initial convec- 
tive motion is in radial direction. However, due to the con- 
tinuity of fluid motion and the non-liner interactions among 
the turbulent eddies, a part of the initial turbulent kinetic 
energy will be converted into the horizontal directions. It is 
known from the theory of turbulence that the correlation 
of turbulent pressure and velocity gradient tends to make 
turbulence isotropic (Rotta 1951). The convection param- 
eter C3 in eq. ((8)l is introduced to describe the tendency of 



restoring the isotropy. The larger ca , the stronger the restor- 
ing force and the turbulence is therefore more isotropic. 
Therefore, the convection parameter ca can be regarded as 
the parameter for the degree of isotropy of turbulence. In 
the convectively unstable zone, the ratio of the radial com- 
ponent of turbulence (wj.) and the horizontal one (wj,) is 
1^1^ « (3 + ca) /2c3 (Deng & Xiong 2006). 

Eqs. (Il])-(|8} are a completed set of dynamical equa- 
tions of stellar envelope structure, which can be used to 
construct the equilibrium envelope model and to calculate 
the radial oscillations. When setting all the terms involv- 
ing time derivatives and velocity Ur in these equations to 
be zero, we will get the equations for equilibrium model. 
When a star does small amplitudes oscillations around the 
equilibrium, all the quantities can be made to be the sum of 
the equilibrium values and the perturbation (oscillations), 
then a Tyler expansion around the equilibrium can be ap- 
plied. By conserving only the linear terms and neglecting all 
terms over the flrst order, we will get the equations of lin- 
ear non-adiabatic oscillations. Details on the equations for 
equilibrium model and linear pulsation, and the method for 
solving these equations are not given here, for those who are 
interested in that, our previous work is referred (Xiong & 
Deng 2001b, Xiong, Deng & Cheng 1998). It needs to be 
noticed that we assumed quasi-isotropy for turbulent con- 
vection in our previous work. In the current, turbulent con- 
vection is more carefully handled, with anisotropy of turbu- 
lence included, i.e. a equation designated to the anisotropy 
of turbulent convection ^ is now added to the equations 
for both the equilibrium model and oscillations. 



3 NUMERICAL RESULTS OF LINEAR 
NON-ADIABATIC OSCILLATIONS FOR 
RED GIANTS 

Making use of the equations introduced in the previ- 
ous section, the linear non-adiabatic oscillations of 10 
sequences of evolutionary models with masses M — 
1, 1.2, 1.4, 1.6, 1.8, 2.0, 2.25, 2.5&2.75M0 are calculated. The 
chemical composition of all these models is solar, X — 0.70 
and Z = 0.02. The evolutionary tracks are made using 
Padova code (Bressan, et al. 1993). For selected stellar 
models along the tracks with given luminosity and effec- 
tive temperature, the corresponding non-local convective 
models of equilibrium envelope are constructed, then lin- 
ear non-adiabatic oscillations for the fundamental through 
11th overtones are further calculated. Special care is taken 
when handling thermodynamical and dynamical coupling 
between convection and oscillation as describe in detail in 
the previous section. Table [T] gives, as an example, the am- 
plitude growth rate r — —2nuji/ujr of linear non-adiabatic 
oscillations for one of the evolutionary model with mass 
AI = I.OM0, where coi and Ur are respectively the imag- 
inary and the real component of the complex frequency of 
oscillations u = icui + uir . 

The equation of state adopted in this work is a simpli- 
fied MHD expression. The main simplification made here is 
that a helium atom is still considered as a quasi-hydrogen 
atom, which makes the calculation for the energy levels of 
helium atomic far easier. For the first ionized helium atoms 
or the high level excitation states of neutral helium atoms. 
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Table 1. Amplitude growth rates for a set of evolutionary models of star with IM© 
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this is no doubt a fairly good treatment. This approxima- 
tion turns out to have large deviation from the real situation 
only for the low excitation states of neutral helium. It is pre- 
dictable that the quasi-hydrogen atom approximation only 
brings into the equation of state and the relevant thermody- 
namic quantities very limited affection, however, it has the 
very good advantage that, instead of interpolated with the 
pre-computed tables, the equation of state and the thermo- 
dynamic quantities can be calculated on the mesh in real 
time. The radiative opacity used in this work is a analytic 
formula to approach the OPAL tables (Rogers & Iglesias 
1992) and the low temperatures opacity tables (Alexander 
et al. 1994). This analytic expression gives a smooth connec- 
tion of OPAL table at intermediate and high temperature 
region and that of Alexander et al. (1994) in low temper- 
atures region. The surface boundary for pulsation calcula- 
tions are set at r = 0.01, and the bottom one is located at 
a temperature of T ~ 5x10^ K. 



3.1 Pulsation stability 

The distribution of the pulsationally stable and unstable 
models in HRD from our calculations is shown in fig. [T] 
where the tiny solid dots are stable models from the funda- 
mental through the 4th overtone (r^ < 0, i = 1, 4); while the 
open circles, crosses, triangles, inverse triangles and squares 
are the unstable models respectively in the fundamental 
mode to the 4th overtone. Two separate instability strips 
can be seen clearly in the plot. The one of them in the mid- 
dle of fig. [l]is the well known S Scuti instability strip which 
is not really the topic of current work. As the 5 Scuti strip 
is calculated using the same scheme, showing it here assures 
the reliability of our theory and the numerical scheme. The 
second one is the pulsating red giant instability strip on the 
upper-right of the figure. All the models between the two 
strips are all pulsationally stable at leeist for modes up to 
the 4th overtone. For RGB stars with low luminosity, the 
fundamental and the first overtone are both stable. These 
stars are unstable only in high order overtones. For high 
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Figure 1. The distribution of pulsationally stable and unstable 
stars in the HRD. The tiny dots are pulsationally stable stars from 
the fundamental up to 4th overtone, the open symbols of circles, 
plus, triangles, inverse triangles and squares are respectively the 
unstable stars in the fundamental up to the 4th overtone. 
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Figure 2. The period— luminosity relation for the evolutionary 
sequence of a IMq star. The horizontal axis is the logarithmic 
period, and the vertical one is the the reddening free Wesenheit 
index (Soszyhski et al. 2004). The circles, plus signs, triangles and 
squares represent respectively the fundamental up to 3rd over- 
tone. 



luminosity and lower effective temperature stars (going up 
along the RGB and AGB), pulsation instability shifts to- 
wards lower overtone and the fundamental mode. 

It follows from table[T]and fig.[T]that all the modes lower 
than the 5th order modes are stable for the M = I.OMq 
RGB stars with logTe > 3.66. This agrees well with the 
observations of Henry et al. (2000). For stars with higher 
masses, the blue edge of the red giant instability strip shifts 
to higher temperature and higher luminosity. For a M = 
3Mq RGB star, the blue edge of instability strip moved to 
logTe ^ 3.68. 

It has to be bared in mind that the location of the 
blue edge of the red giant instability strip depends not only 
on stellar mass, but also on metallicity. For increasing (de- 
creasing) metallicity, the blue edge moves to red (blue) . We 
did similar calculations for intermediate or low luminosity 
stars in globular clusters {Z = 10~*) and we found that 
all RGB stars below the HB are stable in the fundamen- 
tal mode while oscillating in 2-4th overtones, in that case 
the blue edge was shifted to logTe — 3.74 (Xiong & Deng 
2006a,b). The 2 small amplitude short period red variables 
K1098 and VZ1140 found in M15 and M3 (Yao & Qin 1993, 
Yao, Zhang & Qin 1993) are both located near the intersec- 
tion of HB and RGB, and their periods of the fundamental 
mode are both greater than 1 day according to theoretical 
calculations. However, the observed periods are respectively 
0.5155 and 0.35 days, therefore suggesting that they are very 
likely pulsating at 2nd to 4th orders overtone. This agrees 
with the prediction of our work. 



3.2 Pulsation period and constant 

The pulsation period can be derived from pulsation cal- 
culations. It is then straightforward to build a theoretical 
period-luminosity relation. However, to make such relation 
ready to be compared with that of observations is not so 
easy. In fig. (2] the theoretical period-luminosity relation for 
a, M = I.OM0 RGB model is given. The horizontal axis is 
the pulsation period, and the vertical one is the reddening 
free Wesenheit index Wi (Soszyhski et al 2004) 

Wi^I- 1.55 (V - I) 

The temperature to color index conversion and the bolo- 
metric correction is taken from Johnson (1966). It can be 
found that the overall result of our calculations looks very 
similar to the observational one made for the red giant vari- 
ables in Magellanic clouds (Soszyhski et al. 2004). Unfortu- 
nately, quantitatively accurate comparison between our the- 
oretical result and the observational relation is not possible 
for now due to the following reasons, 

(i) The temperature of the RGB stars depends critically 
on the metallicity, low temperature opacity and the mixing 
length parameter a. Larger metallicity and higher opacity, 
or smaller mixing length parameter a, will make the RGB 
move towards lower effective temperature; 

(ii) For low temperature RGB stars, both the effective 
temperature to color index conversion and bolometric cor- 
rection are the worst determined, and this is especially true 
for AGB stars usually having a circumstellar dust envelopes 
which complicated the problem; 

(iii) The stars in the Magellanic clouds are not at all a 
group of stars having the same mass, age and chemical com- 



position. Instead, they have all combinations of all those 
parameters. 

In practice, the ratios between different overtones is a 
far better predictable than the periods, because they are 
less affected by stellar parameters. Figs. [3^ and|3ji are the 
theoretical Petersen diagrams from our calculations. A di- 
rect comparison shows that our results are very similar to 
the observations for the OSARG in LMC (Soszyhski et al. 
2004). Therefore we make a further speculation. There are 
the 4 ridges of period-luminosity relations in their fig. 7. We 
identify that the ridge (or sequence) ai (bi) is the 1st over- 
tone, a2 (b2) is the 2nd overtone, while as (bs) is a mixture 
of the 3rd and 4th overtones, and a4 and b4 are the 5th over- 
tone. In fact, P4 in their paper should be the 5th instead of 
4th overtone, because for the 4th overtone, it is not possible 
to have ratios of P4/P1, P4/P2 and P4/P3 being as small 
as 0.39, 0.56 and 0.76 respectively. However, if one considers 
that P3 in their paper is a mixture of 3rd and 4th overtones, 
while P4 is the 5th one, then a perfect match between the 
observations and our theoretical model. 

Relative to pulsation period, the pulsation constant Q 
changes less, since it is only a slow varying function of stel- 
lar mass, luminosity and effective temperature. Fig.|4^ shows 
the pulsation constants Qi (i=0-3) of the fundamental mode 
through the 3rd overtone as functions of effective tempera- 
ture, calculated for 4 sets of RGB models with masses of 1, 
1.4, 2.0 and 3.0 Mq respectively, and with the same lumi- 
nosity L = 10001/0 . From fig. |4^, the pulsation constants of 
lst-3rd overtones changes relatively smaller than the funda- 
mental mode does. The pulsation constants grow slowly with 
decrease of effective temperature for these overtones. While 
for the fundamental mode, on the contrary. The pulsation 
constants change largely with complex pattern, being a non- 
monotonic function of effective temperature, but showing a 
tendency of slight decrease with increase of mass. 

Fig. |4]d shows the pulsation constant as a function of 
luminosity for 4 sets of models of red giant models with the 
same temperature log Te — 3.55 but the different masses. 
Similarly, the pulsation constants of the lst-3rd overtones 
are almost unchanged, while the fundamental mode shows a 
much larger variation, the pulsation constants increase with 
increase of luminosity of stars. 

To further explore how the pulsation constant changes 
with stellar parameters such as mass, luminosity and ef- 
fective temperature, 24 sets of models of red giant are 
calculated as a supplementary database for this aim. 
They consist of 600 models with diffirent masses (Af = 
1, 1.4, 2.0and3.0A/o), luminosites (logL/L© = 2.6-3.5) and 
effective temperature (logTe = 3.40 — 3.70). we get the fol- 
lowing analytic fitting formula for pulsation constant, 
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m = log M/Mq, 

f = log Te- 3.55, 

Z = logL/LQ -3.0. 

The coefficients in above expressions, (a — /)„, from 
the fundamental to the 3rd overtone are given in table [S] 
Figs. [5K and [5Jd demonstrates the difference A log Qi — 
logQi — logQi.ap (i=0-3) between the calculated pulsation 
constants Qi for the nearly 600 models and their analytic fit- 
ting Qi^ap as functions of effective temperature logT^ ((S^) 
and luminosty ^jp). It is clear from fig.[5]that a pretty good 
fitting is secured in the domain defined by 1 < M /Mq < 3, 
3.4 < logTe < 3.7 and 2.6 < logL/L03.5. 



4 THE EXCITATION MECHANISM FOR THE 
RED GIANT VARIABLES 

By using a work diagram, the excitation and damping on 
the oscillations of stars due to different factors can be clearly 
illustrated and analysed quantitatively. The definition and 
derivation of the integrated work (IW) can be found in the 
Appendix. The total IW Wall is composed of the following 
3 components. 

Wall ^Wp + Wpt + W^.s (22) 

where Wp is the pressure component (including radia- 
tive pressure), Wpt is the turbulent pressure one, while Wvis 
is that of turbulent viscosity, all being functions of depth 
[Mr or log P). In our work, the IW is counted from the stel- 
lar center and integrated outwards, therefore the total IW 
Wall at the surface should be the growth rate of pulsation 
amplitude rj — — 27raJi/cJr., where uji and lj,. are respectively 
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Figure 3. The theoretical Peterson diagram for evolutionary model sequence. Plotted are the period ratios Pj/Pi {i = 0,1,2,3; j = 

i — 1, 5) versus logP^. The 4 panels have different i indices, a) i = the fundamental mode; b) i = 1 the 1st overtone; c) i = 2 the 

2nd overtone; and d) i = 3 the 3rd overtone. The circles, plus signs, triangles and squares are respectively the sequences for M=l, 1.4, 
2.0 and 3.OM0. 



the real imaginary and real part of the complex pulsation 
frequency lu = iuJi + Wr, i.e., 



V = Wa„ (Mo) 



(23) 



When Waii{Mr) (or its one component) increases out- 
wards, this region is therefore a pulsational excitation zone 
(of the corresponding component), otherwise it will be a 
damping zone. 

It is well known that the quasi-Cepheid variables (in- 
cluding Cepheids of populations I and II, RR Lyrae and 
S Scuti stars, this term will be used hereafter) are excited 
by the k mechanism. The existence of the blue edge of the 



instability strip is due to the fact, when stellar surface tem- 
perature is too high, the ionization zones of hydrogen and 
helium that is responsible for k mechanism are very near 
the surface, their mass and heat capacity are too small, and 
they produce far less excitation that is not enough to com- 
pensate the damping at the deep interior, therefore oscilla- 
tions cannot excited for these hot stars. This is shown by the 
plot of IW versus depth in fig. [6]for a hot star (M = 2.0Mq , 
log L/Lq and logTe — 3.9222) located outside the blue edge 
of the 5 Scuti instability strip. It can be seen from fig. [S] 
that convection in the whole envelope is negligible. In the 
deep interiors below the second ionization zone of helium, 
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the overall trend of the integrated work Wp is decreasing 
towards the surface of ster due to the radiation damping. 
Near logP = 8, the Wp curve shows an abrupt ascending 
followed by a quick descending, this is due to the iron ab- 
sorption peak over there. For log P ~ 6.0 — 4.0, Wp increases 
quickly, and it is excited by the k mechanism in the first and 
second ionization zones of helium. There is a little rise up 
at logP ~ 4.0 which is made by the k mechanism due to 



ionization of hydrogen. Nevertheless, all the excitation in 
the ionization zones of hydrogen and helium are not enough 
to compensate the radiation damping arise from the deep 
interiors, this makes Wp < for this model, therefore the 
star is pulsationally stable. 

The ionization zones of hydrogen and helium go deeper 
inside with decrease of effective temperatures of stars, their 
total mass and heat capacity become larger, therefore ex- 
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Figure 6. The integrated work (IW) of the fundamental mode for 
a hot star located outside the blue boundary of the 5 Scuti insta- 
bility strip. The solid line is the integrated work (IW), dashed line 
the opacity, and dotted line the fractional radiation flux Lr/L as 
functions of depth (log P. The 3 heavy horizontal line on the up- 
per left part of the figure represent the locations of the ionization 
zones of hydrogen and the two of helium. 
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Figure 7. The integrated work (IW) of the fundamental mode 
for a warm star within the & Scuti instability strip as a function 
of depth (the solid line). The runs of Wp (the dashed line, gas 
pressure component), Pt (the long dashed line, turbulent pressure 
component), W^is (dotted line, turbulent viscosity component), 
and log Lr/L (dot-dashed line) are also shown. The parameters 
including mass, luminosity and effective temperature of the star 
are labeled. The coupling between convection and oscillations has 
been taken into account. 



citation gets enhanced. When the excitation takes over the 
radiation damping in the deep interior, stars will become 
pulsationally unstable. Fig. [7] depicts IW versus depth for a 
warm star within the 5 Scuti instability strip (M = 2.QMq, 
logL/Lodot = 1.5246 and logT^ = 3.8272). This star pos- 
sesses a relative shaUow convective zones of the hydrogen 
and first helium ionization The second ionization zone of 
helium is radiative in this case. Around log P = 5.5 both 
Wp and Wall increase quickly outwards due to the excit- 
ing of the second ionization zone of helium. The convection 
component (Wpt and Wvis) is much smaller than the exci- 
tation due to radiation. This excitation is also greater than 
the radiation damping in the deep interiors, therefore the 
star become unstable. 

Going to even lower effective temperatures, convection 
zone in stars becomes deeper. When the second ionization 
zone of helium becomes complete convective, the excitation 
due to radiation «;-mechanism reduced largely. In the deep 
interior, however, the thermodynamic coupling between con- 
vection and oscillations works now as a damping for pulsa- 
tion that makes the star stable, this is true at least for the 
fundamental and the first overtone. The IW for a yellow gi- 
ant (M = 2.OM0, logL/L© = 1.3426 and logTe = 3.7368), 
which locates in between the 5 Scuti and the pulsating red 
giant instability trips, is shown in fig. [S] Our calculations 
demonstrate that this star is stable at least in the first 7 low 
order overtones. 

With even lower effective temperature and higher lu- 
minosity, pulsation instability moves to lower order modes. 
In figs. |5^-[S}1 shows IW of the fundamental and the first 
3 overtone modes versus depth for a star on the RGB 
(M = 2.OM0, logL/Lo = 2.3451 and logTe = 3.6301). 
For red giant stars, the pulsation amplitude decreases very 
quickly going downwards from the surface, and tends to 
vanish at the center. The effective pulsation region can be 
seen directly from the IW plot, in which the radiative flux 
is much smaller than the convective one. This infers that 
the gas pressure component of IW, Wp, mainly comes from 
the thermodynamic coupling between convection and oscil- 
lations. In the deep interior (logP > 5), the general trend 
of Wp is to decrease with decreasing log P, meaning that 
the thermodynamic coupling between convection and oscil- 
lations functions as a damping. Near to the stellar surface, 
Wp increases outwards, that is resulted from the radiative 
modulation excitation in the radiative flux gradient zone 
(Xiong, Cheng & Deng 1998). It is clear from figs. 
that, for these red giant stars, the effects of radiation and 
the thermodynamic coupling between convection and oscil- 
lations (Wp) are already much smaller than that of dynam- 
ical coupling {Wpt and Wvis)- For red giants with inter- 
mediate to low luminosity, Wpt (>0) and Wvis (< 0) have 
similar absolute value while having opposite signs. The fun- 
damental and the first 2 overtone modes have Waii < 0, 
therefore are pulsationally stable; for modes higher than the 
3rd overtone, however, Waii > therefore the stars become 
unstable. As shown in table [T] all the low order modes are 
stable for low luminosity red giants which are unstable only 
at the high order overtones. Going to stars of higher lumi- 
nosity and lower effective temperature, pulsation instability 
shifts steadily towards lower order modes, the pulsational 
amplitude growth rates increase very quickly. For the lu- 
minous Mira variables the fundamental mode and the low 
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oder overtones are only linear unstable modes, all the modes 
above 4th order become stable (Xiong, Deng & Cheng 1998) . 
Such a trend is marginally shown in fig. [1] because only the 
first 5 pulsation modes are displayed, and no pulsation sta- 
bility coefficientscan be given. Table [T] is can make up this 
shortcoming, which shows that all modes below 8th over- 
tone are pulsationally stable for all population I stars with 
M = I.OMq and logT > 3.70. It is clearly shown from ta- 
ble[l]that the pulsation instability shifts towards lower order 
modes and amplitude growth rates increase with increase of 
luminosity and decrease of effective temperature of stars. 
Due to numerical difficulties, we unfortunately could not 
construct equilibrium model for high luminosity red giants 
under complete non-local treatment of convection. The most 
luminous and coolest model got by us is No. 40 in the table. 
However, the general trend mentioned above can be seen 
confidently from the table [T] and fig. [T] This trend agrees 
well with our previous work on the luminous long period 
variables (Xiong, Deng & Cheng 1998). The only difference 
between previous and the current works is the constructing 
of the equilibrium model, the non-adiabatic oscillations are 
calculated in exactly the same way: both are treated with 
our non-local convection theory. As discussed in the previous 
section, the amplitude growth rates of SARVs are smaller 
than those of luminous Miras by 1-3 orders of magnitude, 
very special care is required in the current work where all the 
static equilibrium models for the non-adiabatic oscillations 
of SARVs are computed using our complete non-local con- 
vection theory. However, the approximate quasi non-local 
convection models are applied in the previous work for the 
high luminosity LP variables, which is a modification of the 
local convection models. In the convectively unstable zone, 
all the convective quantities in the quasi non-local models 
take the same value as in the local convection models, while 
in the overshooting zone, they are extended from the lo- 



cal models following a analytic asymptote of the properties 
of overshooting outside the unstable zone. Comparing the 
numerical calculations between local and non-local models, 
it shows that within a convectively unstable zone, they are 
rather close to each other, with the sensible difference only 
near the boundaries and in the adjacent overshooting zones 
(Xiong 1989b). For luminous red variables having very large 
pulsation amplitude growth rates, such a approximate treat- 
ment may induce a certain amount of uncertainty to value 
of the growth rates, however, the qualitative results on the 
stability of these stars will not be affected. The numarical 
tests confirm this inference. 

By the way, the so-called frozen convection approxima- 
tion is still favored by some workers (Gong et al. 1995) in 
the study of non-adiabatic oscillations of stars. They assume 
that convection is unchanged during the course of stellar os- 
cillations. Such an approximation is correct only when the 
pulsation timescale is much smaller than that of convective 
motion. Unfortunately, these timescales are usually the same 
order of magnitude within the convective zones. Therefore 
it is clear that this kind of treatment must be incorrect in 
most cases, especially for analyzing the instability of the low 
temperature stars with extended convective envelopes. The 
only chance for the frozen convection approximation to be 
reasonable is for the study of the blue edge of the instability 
strip, because for hot star near the blue edge, convection 
zone is rather shallow so that radiation dominates, and the 
coupling between convection and oscillations has little infiu- 
ence on the stability of these objects. In any case, this ap- 
proximation can explain neither the red edge of the Cepheid 
instability strip, nor the blue edge of the pulsating red gi- 
ant instability strip. We would like to conclude confidently 
that if the coupling between convection and oscillations is 
excluded, namely the frozen convection approximation is ap- 
plied, all stars located to the red side of the blue edge of the 
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Figure 9. Similar to fig. [T] but for a RGB star, a) the fundamental mode; b) the 1st overtone; c). the 2nd overtone; d) the 3rd overtone. 



Cepheid instability strip in the HRD will be pulsationally 
unstable. 

Figs. llOb , and llOb display the IW versus depth for the 
fundamental and the 1st overtone of the same stable star 
as in fig. |8] which is located in between the 5 Scuti and the 
red giant instability strips, but the coupling between convec- 
tion and oscillations is not considered in this calculation. In 
fig. llOb . there are two rapid increase of IW in the surface and 
bottom boundaries of the convection zone. They are resulted 
from the radiative modulation excitations produced in the 
radiative flux gradient regions. Although being a result of 
radiative flux and related to opacity, the radiative modu- 
lation excitation is completely distinct from the radiative 
K mechanism. The former one exists only in the radiative 
flux gradient region whose functionality can be found in our 



previous work (Xiong, Cheng & Deng 1998). Obviously, the 
ionization zones of hydrogen and helium where k mechanism 
operates, are now fully convective, in that the radiative flux 
is only a tiny fraction of the total flux. As shown in fig. llOb ,. 
Wp is virtually a horizontal line in the ionization zones of 
hydrogen and helium. In the deeper interior, radiation ought 
to be a damping, but in the radiative flux gradient region 
just inside the lower boundary of the convective zone, Wp 
shows a abrupt jump. This is a good evidence to show that 
the radiative modulation excitation is not the usual k mech- 
anism. For the 1st overtone (fig. llOb ). Wp shows only a 
single abrupt jump inside the upper boundary of the con- 
vective zone. This is due to the fact that oscillations tend to 
be in further outer layer for the 1st overtone than the funda- 
mental mode, near the lower boundary of convective zone. 
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Figure 10. The integrated work (IW) as a function of depth for the same star as fig. [S] but without talcing into account the coupling 
between convection and oscillations, a) the fundamental mode; b) the 1st overtone. 



the pulsation amplitude is already very small, the radiative 
modulation excitation is therefore much less. 

Figs and [SJd show the IW versus depth for the same 
stellar model as the one described in above paragraph, but 
having the coupling between convection and oscillations 
taken into account. There is no any similarity in the shapes 
of IW shown in fig. |8] and fig. 1101 In fig. [S] the absolute val- 
ues of Wpt and Wvis are both much greater than that of 
Wp, and the Wp curve runs completely in a different form 
from that in fig. 1101 This difference is caused by coupling 
between convection and oscillations. The effects of radiation 
have aready overwhelmed by that of convection for these low 
temperature stars. 



5 CONCLUSIONS AND DISCUSSION 

In this work, we have studied the linear non-adiabatic os- 
cillations of the evolutionary models with M = 1 — 3AIq 
using a non-local time-dependent theory of convection. The 
excitation mechanism for the instabilities are studied. The 
main results can be concluded in the following, 

(i) Two separate pulsational instability strips are found in 
the HRD: the hotter on is the 5 Scuti instability strip, while 
the cooler one is the pulsating red giant instability strip. 
The two strips are intercepted by a region of pulsationally 
stable (at least from the fundamental up to the 4th overtone) 
yellow giants; 

(ii) The 5 Scuti stars and the pulsating red giants (PRGs) 
are two distinct types of variables. The former one is excited 
by the radiative k mechanism and the later one is drived by 
the coupling between convection and oscillations; 

(iii) For the low luminosity PRGs, pulsations are in the 
high order modes. Going to lower effective temperature and 



higher luminosity, pulsationally unstable modes shift grad- 
ually to lower orders from higher orders; 

(iv) In the deep interior of convective zone far away 
from the boundary, generally speaking, the thermodynamic 
coupling between convection and oscillations is a damping 
against pulsation This is the true reason for the existence of 
the red edge of the Cepheid instability strip. Roughly speak- 
ing, when convection zone extend below the second ioniza- 
tion zone of helium, the red edge of Cepheid instability strip 
shows up; 

(v) Obviously, turbulent viscosity is always a damping 
against oscillations, whereas turbulent pressure is a excita- 
tion mechanism. The dynamical coupling between convec- 
tion and oscillations plays a key role in exciting the oscilla- 
tions of PRGs; 

(vi) For the instability of the low temperature stars hav- 
ing extended convective zone, the coupling between convec- 
tion and oscillations is absolutely a dominant factor. When 
the coupling is excluded, neither the red edge of the Cepheid 
instability strip nor the blue edges of the PRGs instability 
strip can be explained, and all stars on the right side of the 
red side of the blue edge of the Cepheid instability strip will 
be pulsationally unstable; 

(vii) The theoretical Peterson diagram given in this work 
is very similar to that of observations by Soszyhski et al. 
(2004). However, based on our calculations of period ratios, 
we would like to guess that very likely the P4 in their work 
is in fact the 5th overtone, while P3 is a combination of the 
3rd and 4th overtones. For the observed modes with period 
ratio close to 1, no explanation can be made at this point. 

PRVs are a group of low temperature red giant stars 
located outside the Cepheid instability strip in the HRD. 
Their exciting mechanism was not known in the past. Based 
on a non-local time-dependent theory of convection, the 
present work provides a natural solution that explains 2 
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types of variables of distinct properties: 5 Scuti stars (and 
all the Cepheid like variables) and the PRGs. In previous 
section, we have addressed the existence of the red edge of 
the Cepheid instability strip, and the presence of another 
PRV instability strip at a lower temperature region in the 
HRD outside the Cepheids instability strip. 

In the General Catalogue of Variable 5'tars(Kholopov 
et al. 1985-1990), the long period variables (LPVs) are clas- 
sified as 3 primary types: Miras, semi-regulars (SRs) and 
iiregulars (L) . Following the statistical study of PRVs in the 
Galactic field and stellar group (stellar clusters and asso- 
ciations), Eggen (1973a,b, 1975, 1977) made another clas- 
sification according to pulsation amplitude, which is again 
in 3 major types: Large-amplitude red variables (LARVs: 
/S.V > 2mag.), medium-amplitude red variables (MARVs: 
0.3 < AV < 2mag.) and small amplitude red variables 
(SARV: AV < 0.5mag.), where the ones with visual am- 
plitude less than 0.2 mag. and very short period (20- 
40days) may form a individual group. Based on the MACHO 
database. Wood (2000) discovered that the PRGs in LMC 
displayed 5 ridges (or sequences) of A, B, C, D and E in the 
K — log P diagram. The Miras, located in sequence C, have 
the largest amplitude and for sure represent the fundamen- 
tal mode. Stars falling into sequence B are SRs possessing 
smaller amplitudes. In sequence C, one finds the stars with 
relatively smaller amplitudes {AV < 1.5mag), also being 
SRs. Stars with very small amplitudes {AV < 0.2) falling 
in sequence A are very likely the SARVs as classified by 
Eggen. Sequences B and A may correspond to the lst-3rd 
overtones (Wood, 2000). Soszyiiski et al. (2004, 2005) con- 
firmed the discovery of Wood using OGLE-II and OGLE- 
III data, and they further pointed out that the OSARGs 
(OGLE small amplitude red giants=SARVs) and the LPVs 
(Miras) are two different types of pulsating stars. On the 
Period-Luminosty (PL) diagram (periods-W/), the primary 
periods of OSARGs fall in sequences A and B, and none 
of the second and the third dominant period-W/ relation 
of OSARGs fall in sequences C; On the contrary, the pri- 
mary period-W/ relation fall in sequences C and B, and 
none of the second and the third dominant period-W/ re- 
lation goes in sequence A. Based on classification as such, 
in the PL diagram of Wood's ridges A and B of OSARGs, 
can be divided into 4 separate ridges ai, a,2, aa and sl4 for 
OSARGs brighter than the tip of RGB, and bi, b2, bs and 
b4 for OSARGs fainter than the tip of RGB. According to 
the theoretical Peterson diagram, it is confident to believe 
that ridge C is the fundamental model, ai (bi) and a,2 (b^) 
are respectively the 1st and the 2nd overtones. Whereas, as 
and a4 are not likely to be single mode sequence: as (bs) 
is a mixture of the 3rd and the 4th overtones, and a4 (b4) 
is a mixture of the 4th and the 5th overtones. Otherwise, 
as discussed in section 3, it is difficult to understand the 
observational Peterson diagram (Soszyhski et al. 2004). As 
shown in fig. O the period ratios are hard affected by stellar 
mass, following similar arguments, they are neither affected 
by metallicity. Because metallicity determines only the lo- 
cation of the instability strip and the ML relation. As both 
mass and metallicity hardly infiuence the period ratios, our 
results for mode identifications should be reliable. 

Observational evidence provided by Wood and 
Soszynski et al. demonstrate that there is a tight link be- 
tween the properties of PRGs and their oscillation modes. It 



is clear that primary period of Miras are in the fundamental 
mode (Soszyhski et al. 2004), the primary period of SRs are 
either in the first overtone or the fundamental mode, the sec- 
ond the third dominant periods of them should not be higher 
than the 4th overtone. On the contrary, most of the of the 
SARGs are pulsating in the 3rd-4th overtones (~ 50%) or 
in the 2nd overtone (~ 20%), and their fundamental mode 
is pulsationally stable. How to inspect theoretically such ob- 
servations? 

Based on the theory of oscillations, when neglecting 
both rotation and magnetic field, the pulsation properties of 
stars depend primarily on luminosity and effective tempera- 
ture, and secondarily on mass and metallicity. It follows from 
table[T]that, for the faint RGs, pulsation instability starts in 
high order overtones, and it shifts to lower orders with in- 
crease of luminosity and decrease of effective temperature, 
while the pulsation amplitudes also grows. By comparing the 
pulsational instability for different stellar masses, as shown 
in fig.[T] going towards higher masses, the blue edge of the in- 
stability shifts to higher luminosity and higher temperature. 
This theoretical results agree well the observations. Eggen 
(1977) studied the distribution in Mi,o;-(R-I) diagram for 
PRVs in the Galactic field and stellar groups, and found 
out that the SARVs are indeed located in the high temper- 
ature (0.8 < {R — I)o < 1.6) and low luminosity region, 
the LARVs stay in the low temperature and high luminos- 
ity area, whereas MARVs stay in the intermediate region. 
Henry et al. (2000) also shown that all the K5-M0 giants are 
SARGs, they are located in the high temperature and low lu- 
minosity region in the HRD. Percy & Parkes (1998) studied 
the pulsation modes for 13 SARGs with well-determined pe- 
riod and accurate Hipparcos parallaxes, and claimed that 9 
of the stars are pulsating in the 1st or the 2nd overtone, 1 in 
the 3rd overtone, while only 3 in the fundamental or the 1st 
overtone. Soszyhski et al.'s observational study shown that 
non of all their OSARGs is oscillating in the fundamental 
mode. Therefore, it can be conclude that our theoretical pre- 
dictions agree pretty well the observational trend. However, 
it is difficult to understand that, on the log P-M/ diagram 
of Soszyhski et al. (2004), the OSARGs extends almost all 
the way to the brightest end. We have not yet understant 
OSARGs. Why they are so different from the Mira and SR 
variables at the same luminosity? It can be seen from fig. [T] 
that the blue edge of PRGs instability stip move toward 
higher temperature and luminosity with increase of stellar 
mass. OSARGs, we suppose, are more massive and hotter 
than Mira and SR variables with the same luminosity. The 
brightest OSARGs may be a group of younger massive red 
giants. If this guess is true, there shoud not be lumilous 
OSARGs in the old clusters. 
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where Pt = px^ is the turbulent pressure, uj is the com- 
plex frequency of oscillations defined as, 

U! = LUr + iuJi, (A2) 

putting HA2|) into (|A1|) . solving the equation for u, we 



have, 
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The next step is multiplying eq. (|A3|) 
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with 



TT (Sr* / r) dMr , then integrating it for dMr from the 
center to the surface. Baring in mind that the relative 
pulsational amplitude at the center is zero and the pressure 
and density vanish at the surface, and Separating the real 
and imaginary parts, we have, 
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where Im means to take the imaginary of the term in 
the brackets. —2Tiuji/ur is the amplitude growth rate of os- 
cillation within a cycle ry. 
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APPENDIX A: THE INTEGRATED WORK 

Using the eq.[2]in the text, its linear pulsation equations can 
be give as. 



As from eq. (|A5|) . r/ can be divided into a gas pressure 
component and a Renold's stress component. We can define 
the IWs as functions of depth, Wp{Mr) and WdyniMr) , to 
describe the contributions of the two components, 
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The gas component of the IW, Wp(Mr), contains both 
the contribution of radiative energy transfer and the con- 
vective one. Starting from the energy equation in the text, 
eq. [31 we can have, 
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where Lr, and Lt are the luminosities corresponding 
respectively to the radiative, the convective and turbulent 
energy fluxes, ei is the nuclear energy generation rate (with 
that of neutrino subtracted). Due to the fact that oscillations 
drops off abruptly towards stellar center for giant stars, the 
calculation of oscillating stellar model can be simplified by 
doing the job only for the models of envelope. Under such a 
condition, the contribution of energy generation rate can be 
safely neglected. Putting eq. (|A9P into eq. <\A7\i . the contri- 
bution of radiative energy transfer can be singled out from 
that of convective energy transfer, 



Wp (Mr) = Wpr (Mr) + Wth,™ {Mr 
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Re means taking the real part of the terms in the brack- 



ets. 



Wpr is the contribution of radiative energy transfer, 
while Wtherm IS that of convective energy transfer, i.e. the 
thermodynamic coupling between convection and oscilla- 
tions. Wdyn is that due to Renold's stress, the dynamic cou- 
phng between convection and oscillations. The second and 
the third terms in the dynamic equations ([5]) and ((Sjl for 
turbulent Renold stress in the text, represent the turbulent 
viscosity. Solving for SPt and 511^'^ from the linear pulsation 
equations of eqs. ([5} and ((8)l , it is not difficult to isolate the 
viscosity term which is then put in eq. (|A8|I . and we can 
have the turbulence viscosity component of the IW, 
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where Tc is the dissipation time scale of turbulence, 
which is also the inertia time scale of turbulence, expressed 
as, 
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Subtracting eq. (|A8|) with eq. (IA13|) . one has the tur- 
bulent pressure component of the IW, Wpt, 



Wpt {Mr) = Wdyu {Mr) - W r^s {Mr) . 



(A16) 



Summing up the four components, Wpr, Wtherm, Wpt 
and Wvis, we have the total IW, 



Wall (Mr) = Wpr {Mr) + 

Wtherm {Mr) + Wpt {Mr) + Wvis {Mr) , 



(A17) 



On the right hand side of eq. (|A17|I , only the first term 
is the contribution of radiative energy transfer, all the other 
3 terms are due to convection, where Wtherm is the contri- 
bution of the thermodynamic coupling between convection 
and oscillations, Wpt is that due to turbulent pressure, and 
Wvis is that of turbulent viscosity. Under such definitions, 
the linear amplitude growth rate of oscillations, rj, will equal 
the surface value of Waii {Mr), 



7? = -2TTljJi/uJr = Wall {Mq) . 



(A18) 



It is clear from eq. (|A13|I that the turbulent viscosity 
component Wvis of the IW is always negative, therefore it 
is a dumping against oscillations. While the effects of the 
turbulent pressure component Wpt and turbulent thermal 
convection component Wtherm depend very much on the 
phase of the variations of turbulent pressure Pt and turbu- 
lent thermal convection Lc = 4nr^ pcpTV relative to that 
of stellar oscillations {Sp/p, or 5r/r). It follows from the 
dynamic equations of turbulent convection (eqs. [T}{S} that, 
in general, the relative variations of convection quantities 
6x/x, SZ/Z, 6V/V and 5n"/n" always lag behind density 
variation 5p/p (or -Sr/r) inside a convective zone. Under 
local convection approximations, this argument was proved 
analytically in our previous work (Xiong 1977). More general 
oscillation calculations using non-local convection treatment 
further proved this statement (Xiong, Cheng & Deng 1998, 
Xiong, Deng & Cheng 1998). Therefore one can conclude 
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from eqs. ^Mi & IA12|I that Wpt > and Wtherm < gen- 
eral hold inside a convective zone, i.e. turbulent pressure 
is a excitation mechanism for oscillations, while turbulent 
thermal convection (the thermodynamic coupling between 
convection and oscillations) is a damping one. Such a theo- 
retical analysis agrees perfectly with what have been shown 
in figures [IHS] for Wpt, Wtherm and Wvia- 



